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The conjecture that several recently observed mesons have a structure, which is not dominated by 
an ordinary quark-antiquark pair, but by a four-quark structure, is being actively investigated both 
theoretical and experimentally. Such a state may be characterized as a mesonic molecule or as 
a diquark-antidiquark pair. Lattice QCD provides a theoretically sound framework to study such 
states. To quantitatively investigate the internal structure of such mesons, one needs to precisely 
compute correlation matrices containing several interpolating operators including two and four 
quarks. Here we discuss certain technical aspects of such correlation matrices suited to study 
tetraquark candidates with J p = 0 + and flavor structure qxqiq^q-i, e.g. the «o(980) meson, the 
D* 0 meson and some of the charged ccX states. Some numerical results for the «o(980) meson 
are presented. 
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1. Motivation 

Our understanding of the light scalar meson sector (quantum numbers J p = 0 + ) is still incom¬ 
plete [1, 2, 3]. The observed mass ordering of a, K, /o(980) and no(980) is inverted compared to 
expectations from conventional qq quark model calculations. Moreover, when assuming a qq struc¬ 
ture, it is hard to explain the mass degeneracy of /o(980) and no(980). However, the interpretation 
of this meson sector assuming a four-quark structure is rather straightforward. The expected mass 
ordering is consistent with experimental results and due to an additional ss pair the degeneracy of 
/o(980) and «o(980j is not surprising (cf. Figure 1). 

Several lattice QCD studies of the light scalar mesons have been published in the last couple of 
years [4, 5, 6, 7, 8, 9, 10]. In this work we continue our investigation of the ao(980) meson [11, 12, 
13, 14, 15, 16, 17] and briefly discuss the applicability of our methods and codes to other systems, 
the D* 0 meson and some of the charged cc X states. In particular the investigation of tetraquark 
candidates with two heavy and two light quarks seems promising, since recent computations in the 


static limit predict a bound state [18, 19]. 

2. Interpolating operators and the correlation matrix 

Our investigations are based on a 6 x 6 correlation matrix 

Cjk{t) = (g^ {t2)G k \t\)^ , t = t2—t\. ( 2 . 1 ) 

The interpolating operators G* generate quantum numbers I(J P ) = 1(0 + ), 

0 1 = ^=L(£k x ) (2-2) 

^2 = fy KK, point = £ (2.3) 

^3 = fir,,*, point = £ (V x75 , x ) (J x/5Mx ) (2.4) 

G 4 = G Qa =^eabc{sxjb(Cy^)^^e a d e (i^ d (Cy^)s x Jj (2.5) 

^5 = ff KK, 2-part = £ (d y ^ y ) (2.6) 

x ?y 

= 0**' ^ =£ (Wx) (dy75«y) , (2-7) 

x ?y 


where C is the charge conjugation matrix. The operator G q ^ generates a standard quark-antiquark 
state, while all other operators generate four-quark states. G KK ’ pomt and pomt are of mesonic 
molecule structure (KK and r\ s n), while corresponds to a diquark-antidiquark pair (we use 
the lightest (anti)diquarks with spin structure Cf 5 [20, 21, 22]). These three operators are intended 
to model the expected structures of possibly existing four-quark bound states, i.e. of tetraquarks. 
The remaining two operators G KK - 2 " part and G^ 11 - 2 " part independently generate two mesons (K + K 
and 77 v + 71) and, hence, should be suited to resolve low-lying two-particle scattering states. 

In Figure 2 the correlation matrix Cjk is shown in a graphical way in terms of diagrams, where 
each line represents a quark propagator. Clearly there are disconnected diagrams and diagrams, 


2 



Computation of correlation matrices for tetraquark candidates with J p = 0 + 


Joshua Berlin 


mass mass mass 


1 = 0,1 


o o 

1 = 0,1 

ss 

/ 

= 0,1 

_ Issl _ 

1 = 1/2 

- - 

K 

1 = 1/2 

Is si 

I 

= 1/2 

Ills slU 

1 = 0 

— 

O 

To 

II 

o 

_ Jj_ _ 

r 7? 

I 

= 0 

UU 


-1 -1/2 0 +1/2 +1 



-1 -1/2 0 +1/2 +1 



-1 -1/2 0 +1/2 +1 


Figure 1: Spectrum of light scalar mesons ( J F = 0 + ). Experimental results (left); Theoretical expectation 
within the conventional quark model (qq structure) (center); Theoretical expectation assuming a four-quark 
structure (right). 

where quarks propagate within a timeslice. Computing these diagrams precisely is very challenging 
(cf. section 3 and [17]). 
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Figure 2: Diagrammatic representation of the 6x6 correlation matrix Cjk containing the interpolating 
operators (2.2) to (2.7). The diagrams only reflect the spatial quark propagation. Even though elements with 
j,k £ {2,3,4} appear to be identical, they differ in their spin and color structures. 

Note that the correlation matrix given in Eq. (2.1) is also suited to studying the D* 0 meson, 
when replacing the quark flavors according to u —> s, d —> c and s -+ u.d. Similarly, one can explore 
certain charged charmed tetraquarks, when replacing s c. 

3. Techniques to compute the correlation matrix elements 

A number of diagrams of the correlation matrix shown in Figure 2 are challenging to com¬ 
pute, in particular disconnected diagrams or diagrams, where quarks propagate within a times¬ 
lice. Our strategy is to combine several standard techniques including (a) fixed-source propagators, 
(b) stochastic timeslice-to-all propagators, (c) the one-end trick and (d) sequential propagators (cf. 
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[17] and references therein). For most diagrams there are a couple of possibilities, where in many 
cases the most efficient combination with respect to the required HPC resources and the obtained 
statistical accuracy can only be identified after extensive numerical tests and/or analytical consid¬ 
erations. Examples dealing with C\\ and the disconnected diagram of 645 can be found in [17]. 
Below we briefly discuss another example, C 15 . A complete and detailed discussion of the full 
correlation matrix will be part of an upcoming publication. 

The matrix element 

Ci 5 (0 = Yj trrG i (z,fi;y,fi) t y 5 G“(x,t 2 ;z,ti) t 75G :d (x,t2;y,ti))') (3-1) 

(quark propagators are denoted by G q ) can e.g. be computed by combining the one-end trick and 
the technique of sequential inversions. The one-end trick realizes an implicit summation over a 
timeslice, which can be either at t\ (case (a)) or at ti (case (b)), 

(a) : -(I x , z tr^[0 i (z)5- O 4 1 ] t 75G"(x,t2;zTi) T i 75 f / W)) -> -(L tr (v / ’“ A W + 75^W) ) (3.2) 

seq. inversion 

(b) : -(Ly,ztr(</> d (7) t 75p"(y,ti;zTi)750"(z))) -+ -(^tr(</» <i (y) + 75V /j/ “(7))) (3-3) 

seq. inversion ' 

(cf. also Figure 3). At first glance both options seem to be equivalent. Note, however, that the 
computation of the sequential propagator (the blue lines in Figure 3), requires for case (a) only 
2 inversions of the Dirac matrix, while for case (b) t max /a + 1 inversions are necessary to obtain 
Ci 5 (t) for 1 < t < f max . Clearly, proceeding as in (a) is more efficient than proceeding as in (b). 





Figure 3: Two possibilities to compute C 15 using the one-end trick and sequential propagators. Both meth¬ 
ods show comparable errors, but method (b) requires with t/a+ 1 inversions a more computational effort. 


4. Numerical results 

At the moment we are performing computations as described in section 2, section 3 and 
Ref. [17] using around 500 Wilson clover gauge link configurations with 2+1 dynamical quark 
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Figure 4: (top left) Effective masses from the 2 x 2 correlation matrix containing fy KK -v ,mn and 4 , ' ,v7r ' p0inl . 
(top right) Effective masses from the 4x4 correlation matrix containing additionally fj KK - 2 -P' dn anc | 
2-part . ( r |’| le horizontal lines correspond to expected energies of the two-meson states K 4 K and t] s + n, 
where both mesons are at rest and with one relative quantum of momentum.) (center left and right) Squared 
eigenvector components obtained from the 4x4 correlation matrix, (bottom) Effective masses from the 
5x5 correlation matrix containing additionally 


flavors generated by the PACS-CS collaboration [23]. The lattice size is 32 3 x 64, the lattice spac¬ 
ing a ~ 0.09fm and the light u/d quark mass corresponds to m n ~ 300MeV. Since computations 
are ongoing and statistical accuracy is currently rather limited, we neglect in a first analysis dia¬ 
grams of the correlation matrix, where quarks propagate within a timeslice, e.g. the diagram shown 
in Figure 3 and all diagrams with closed quark loops. Consequently, C\j = Cj \ = 0 for j = 2,... ,6, 
i.e. the correlation matrix shown in Figure 2 is reduced to a 5 x 5 matrix containing the interpolating 
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operators 0 1 ,..., . 

We determine effective masses Ef ( and corresponding eigenvectors v„ by solving the standard 
generalized eigenvalue problem 


C(t)\„(t,t r ) = h,{t,t r )C(t r )v n (t,t r ) 


t large „ P ff , , 1 , 

E n = E n n (t,t r ) = -In 


( hn(tfr) \ 
\ hi(f T Ah) / 


(4.1) 


with t r = a. 

Effective masses from the 2x2 correlation matrix containing only the mesonic molecule op¬ 
erators ff KK <v° m and A ' 1 ' 71 ’ 110 " 11 are shown in Figure 4 (top left). There are two states at around 
lGeV, i.e. roughly consistent with the expected mass of the ao(980) meson [24] as well as with 
the two-meson states K + K and T| v + n in our lattice setup. 

Similarly, effective masses from the 4x4 correlation matrix containing additionally the two- 
meson operators (f KK d-v®u an d ffhs^d-part arc s ( lown j n Figure 4 (top right). The same two low-lying 
states around 1 GeV are obtained, but the corresponding plateaus are reached at smaller temporal 
separations t and exhibit less statistical fluctuations. The energies of the second and third excita¬ 
tions are significantly larger, consistent with the expected energies of two-meson states with one 
relative quantum of momentum. This indicates that the two low-lying states are of two-meson 
type, i.e. neither of the two corresponds to a bound four-quark state, which can be identified with 
the «o(980) meson. This is further supported by the squared eigenvector components of Vo and vi 
obtained from the 4 x 4 correlation matrix, which are also shown in Figure 4 (center left and right). 

Results of an analogous analysis of a 5 x 5 correlation matrix containing also the diquark- 
antidiquark operator are shown in Figure 4 (bottom). Since no additional state at around 
1 GeV appears, the conclusions are essentially the same as for the 4 x 4 correlation matrix: It thus 
seems that when neglecting quark propagation within a timeslice, and given the unphysical value 
of the pion mass we are investigating, the no (980) meson is not a rather stable four-quark state 
of mesonic molecule or diquark-antidiquark type. (Cf. also [12] for a more detailed very similar 
discussion.) 


Acknowledgments 

M.W. and J.B. acknowledge support by the Emmy Noether Programme of the DFG (German 
Research Foundation), grant WA 3000/1-1. The work of M.G. was supported by the European 
Commission, European Social Fund and Calabria Region, that disclaim any liability for the use 
that can be done of the information provided in this paper. 

This work was supported in part by the Helmholtz International Center for FAIR within the 
framework of the FOEWE program launched by the State of Hesse. 

Computations have been performed using the Chroma software library [25]. Calculations on 
the FOEWE-CSC high-performance computer of Johann Wolfgang Goethe-University Frankfurt 
am Main were conducted for this research. We would like to thank HPC-Hessen, funded by the 
State Ministry of Higher Education, Research and the Arts, for programming advice. 


6 




Computation of correlation matrices for tetraquark candidates with J p = 0 + 


Joshua Berlin 


References 

[1] J. R. Pelaez, PoS ConfinementX (2012) 019 [arXiv: 1301.4431 [hep-ph]]. 

[2] J. R. Pelaez, AIP Conf. Proc. 1606, 189 (2014). 

[3] C. Amsler et al., “Note on scalar mesons below 2GeV” (2014) 
[http://pdg.lbl.gov/2014/reviews/rpp2014-rev-scalar-mesons.pdf]. 

[4] C. Bernard, C. E. DeTar, Z. Fu and S. Prelovsek, Phys. Rev. D 76, 094504 (2007) [arXiv:0707.2402 
[hep-lat]]. 

[5] C. Gattringer et al., Phys. Rev. D 78, 034501 (2008) [arXiv:0802.2020 [hep-lat]]. 

[6] S. Prelovsek, AIP Conf. Proc. 1030, 311 (2008) [arXiv:0804.2549 [hep-lat]]. 

[7] K. F. Fiu, AIP Conf. Proc. 1030, 305 (2008) [arXiv:0805.3364 [hep-lat]]. 

[8] M. Wakayama and C. Nonaka, PoS FATTICE 2012, 276 (2012) [arXiv: 1211.2072 [hep-lat]]. 

[9] S. Prelovsek, F. Feskovec, C. B. Fang and D. Mohler, Phys. Rev. D 88, 054508 (2013) 

[arXiv: 1307.0736 [hep-lat]]. 

[10] M. Wakayama et al., Phys. Rev. D 91, 094508 (2015) [arXiv: 1412.3909 [hep-lat]]. 

[11] J. O. Daldrop et al. [ETM Collaboration], PoS FATTICE 2012, 161 (2012) [arXiv: 1211.5002 
[hep-lat]]. 

[12] C. Alexandrou et al. [ETM Collaboration], JHEP 1304, 137 (2013) [arXiv: 1212.1418]. 

[13] M. Wagner et al. [ETM Collaboration], PoS ConfinementX, 108 (2012) [arXiv: 1212.1648 [hep-lat]]. 

[14] M. Wagner et al. [ETM Collaboration], Acta Phys. Polon. Supp. 6, no. 3, 847 (2013) 

[arXiv: 1302.3389 [hep-lat]]. 

[15] M. Wagner et al., PoS LATTICE 2013, 162 (2013) [arXiv: 1309.0850 [hep-lat]]. 

[16] M. Wagner et al., J. Phys. Conf. Ser. 503, 012031 (2014) [arXiv: 1310.6905 [hep-lat]]. 

[17] J. Berlin et al., PoS LATTICE 2014, 104 (2014) [arXiv: 1410.8757 [hep-lat]]. 

[18] P. Bicudo, K. Cichy, A. Peters, B. Wagenbach and M. Wagner, Phys. Rev. D 92, no. 1, 014507 (2015) 
[arXiv:1505.00613 [hep-lat]]. 

[19] A. Peters, P. Bicudo, K. Cichy, B. Wagenbach and M. Wagner, arXiv: 1508.00343 [hep-lat]. 

[20] R. L. Jaffe, Phys. Rept. 409, 1 (2005) [hep-ph/0409065], 

[21] C. Alexandrou, P. de Forcrand and B. Lucini, Phys. Rev. Lett. 97, 222002 (2006) [hep-lat/0609004], 

[22] M. Wagner and C. Wiese [ETM Collaboration], JHEP 1107, 016 (2011) [arXiv: 1104.4921 [hep-lat]]. 

[23] S. Aoki etal. [PACS-CS Collaboration], Phys. Rev. D 79, 034503 (2009) [arXiv:0807.1661 [hep-lat]]. 

[24] K. A. Olive et al. [Particle Data Group], Chin. Phys. C, 38, 090001 (2014). 

[25] R. G. Edwards et al. [SciDAC, LHPC and UKQCD Collaborations], Nucl. Phys. Proc. Suppl. 140, 

832 (2005) [hep-lat/0409003]. 


7 



